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I.  Introduction 


Producing  a  plot  of  collected  probe  current  versus  an  applied  potential  sweep  (the  probe 
characteristic)  is  a  standard  plasma  diagnostic  technique  for  determining  electron  density 
and  temperature  in  plasmas.  After  the  original  work  in  this  area1,  there  have  been 
numerous  fundamental  improvements  and  inferences  from  those  basic  ideas.2’3 
Effectively,  the  well-known  method  relies  on  applying  a  very  low  frequency  signal  to  a 
probe  so  that  both  electrons  and  ions  can  fully  respond  to  changes  in  the  driving  voltage.4 
The  dc  characteristic,  Io(Vo),  depends  on  probe  geometry  and  is  used  to  find  the  plasma 
density  via  a  curve  fit,  where  Iq  is  the  dc  probe  current  and  Vo  the  applied  voltage  with 
respect  to  the  floating  potential.  With  these  definitions,  the  dc  resistance  is  given  by  Rjc 
=  Vo/  Io.  The  fitting  procedure,  either  to  the  electron  or  the  ion  saturation  portion  of  the 
characteristic,  has  accompanying  uncertainties  both  with  respect  to  probe  geometry,  and 
errors  associated  with  the  fitting  procedure  itself.  This  can  render  fitting  techniques 
questionable,  particularly  when  probe  surfaces  become  contaminated5  or  when  there  are 
two-temperature  distributions.  Another  issue  associated  with  curve  fitting  is  that  the 
electron  and  ion  saturation  currents,  when  used  to  determine  electron  density,  often 
produce  inconsistent  results.  ’  ’  This  is  an  issue  yet  to  be  resolved  in  practical  probe  data 
analysis  and  is  also  the  subject  of  current  investigations  at  NRL. 

Another  time-honored  method  of  measuring  electron  density  is  based  on  the  rf  impedance 
probe.  This  technique  is  based  on  applying  an  ac  signal  which  is  swept  up  to  and  beyond 
the  plasma  frequency.  A  strong  resonance  occurs  at  the  plasma  frequency  and,  from  that 
resonance,  the  plasma  electron  density  is  readily  determined.  The  applicability  of  this 
technique  to  measurements  from  space  vehicles,  where  motion  of  the  vehicle  can  be 
comparable  to  thermal  ion  speeds,  was  noticed  early  on9  and  has  been  employed  over  the 
years  in  a  variety  of  ionospheric  measurements.  10,11,12  Temperature  is  often  determined 
in  these  cases  by  using  an  accompanying  Langmuir  probe.  ~  With  these  two 
measurements  the  Debye  length  can  be  obtained  and  the  sheath  size  estimated. 

The  rf  impedance  probe  technique  developed  at  the  Naval  Research  Laboratory  is  useful 
for  a  number  of  reasons13,14;  one  of  the  most  important  is  the  unique  detennination  of  the 
bulk  plasma  frequency,  cope,  and  therefore  the  electron  density.  The  network  analyzer 
used  in  these  measurments  can  be  calibrated  to  eliminate  cabling  impedances  and  it 
returns  both  the  real  and  imaginary  parts  of  the  complex  impedance  from  the  reflection 
coefficient.  At  resonance,  the  imaginary  part  of  the  impedance  goes  to  zero  while  the 
real  part  reaches  a  peak,  indicating  maximum  transfer  of  power  to  the  load.  These  two 
conditions  serve  as  a  dual  basis  for  detennining  cape.  The  rf  measurements  are  far  less 
sensitive  to  the  effects  of  surface  contamination  and  secondary  electron  emission  than  the 
conventional  voltage  sweep  analysis. 

As  another  important  application  of  the  network  analyzer  technique,  the  most  recent 
efforts14  were  concerned  with  demonstrating  the  existence  of  collisionless  resistance  (CR) 
in  the  sheath  of  a  spherical  probe  when  driven  by  a  small  rf  signal.  The  CR  was  shown  to 
be  inversely  proportional  to  the  plasma  density  gradient  at  the  location  where  the  applied 
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angular  frequency  equals  the  electron  plasma  frequency,  cope .  The  expression  for  the  CR 
derives  from  a  singular  integral  for  the  complex  plasma  impedance  and  is  given  by, 


R(a>)  = 


n  r  da>pe, 

2e0a}  dr  J'='b 


(1) 


and  is  valid  for  cape(ro)  <  co  <  copb  where  <x>pe(r0)  is  the  plasma  frequency  at  the  probe 
surface,  r0,  and  copb  is  the  bulk  plasma  density.  From  this  basic  expression  we  are  able  to 
construct  the  sheath  density  profile  by  integrating  Eq  (1)  over  frequency  and  inverting. 
Further  details  of  that  work  and  the  origin  of  Eq.  (1)  are  found  elsewhere.14,15 

The  frequency  range  in  this  work  is  applied  such  that  copi«  co  <cope(r0)  where  copi  is  the 
unperturbed  ion  plasma  frequency.  The  ions  are  unable  to  respond  to  changes  in  ac 
voltage  but  the  electrons  are.  Then,  unlike  the  case  where  a  very  low  frequency  is 
applied  (co  <  copi  ,  cope  ),  the  ac  current  will  have  no  contribution  from  the  ions.  This 
implies  that  the  ac  resistance  is  only  a  function  of  the  change  in  electron  current  with 
applied  ac  voltage.  In  this  region  we  are  able  to  derive  expressions  for  the  electron 
temperature,  the  sheath  capacitance  and  the  plasma  potential.  Without  an  ion  response 
we  are  effectively  including  only  the  electron  contribution  to  the  I-V  characteristic. 

We  demonstrate  the  application  of  this  experimental  technique  for  three  small  spherical 
probes  for  differing  environments  and  plasma  conditions.  In  Section  II  we  describe  the 
theory  and  assumptions  associated  with  the  calculations.  In  Section  III  we  outline  the 
experimental  arrangement  in  the  laboratory  and  describe  the  experimental  techniques 
and,  in  Section  IV  we  compare  the  values  of  electron  temperature  detennined  using  this 
method  with  measurements  made  by  a  conventional  Langmuir  probe  analysis  for  each  of 
the  spheres. 

II.  Theory 

a.  Rac  and  Te  determination 

When  we  apply  a  small  ac  voltage  to  a  probe  relative  to  the  dc  bias,  Vac«Vo,  additional 
current  will  flow  to  the  probe.  The  ac  impedance  associated  with  this  addition  is  then 
Zac=  Vac  /  he  and  it  becomes  purely  resistive  as  co  — >  0.  If  we  let  total  probe  current  and 
voltage  be  given  by  V(t)  =  Vo+  VaCe!COt  and  I(t)=Io  +  IaCecot  and  assume,  once  again,  that 
V(t)-Vo  and  I(t)-  In  are  small,  it  follows  that  Zac  =  [V(t)-VfJ  /  [  I(t)-  In]  ~  dV/dl,  where 
I(V)  is  the  dc  probe  characteristic.  The  frequency  region  of  interest  for  this  work  involves 
the  range  where  the  ions  are  not  able  to  respond  to  the  applied  signal  (co»ooPi);  in  that 
case  we  eliminate  ion  current  contributions  and,  in  the  absence  of  displacement  current, 
total  current  is  carried  by  the  electrons  only  or  Zac  ~  dV/dIe.  By  avoiding  any  rf  ion 
current  contribution,  we  are  indeed  observing  dV/dIe  as  opposed  to  dV/dl.  In  addition, 
we  wish  to  avoid  resonance  effects  predicted  by  the  earlier  collisionless  theory  and  so 
require  that  co  <  cope(r0 )  where  r(l  is  the  probe  radius.  As  already  stated,  this  condition 
assures  that  we  are  investigating  the  region  in  frequency  space  where  there  is  no 
contribution  from  the  collisionless  resonance  (/.£.,  there  is  a  sharp  cutoff  since  the  CR 
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occurs  only  when  the  applied  frequency  is  equal  to  the  local  plasma  frequency;  cope(r0 ) 
being  the  minimum  plasma  frequency.  According  to  the  electron  and  ion  momentum 
equations,  neglecting  the  ac  ion  current  is  justified  to  order  nme/neM,  where  me  and  M  are 
electron  and  ion  masses  and  n,  is  ion  density.  To  a  good  approximation  the  electron 
density  obeys  the  Boltzmann  relationship,  ne(r)=no  ep{r)/Te  where  <p(r )  <  0  is  the  local 
potential  with  respect  to  the  unperturbed  plasma  and  Te  is  the  electron  temperature.  Thus, 
although  ne  «  n,  near  the  probe,  copi«  cope  as  long  as  0  <  |  V(/Te  \  <  9  where  Vo  is  here 
defined  as  the  dc  bias  with  respect  to  the  plasma.  The  ac  resistance  is  then  given  in 
Gaussian  units  by, 


Rac=(die/dv0yl 


r V2 

V  ro  J 


-eVn  I  T„ 


(2) 


provided  copi«  co  <  cope(ro).  In  Eq.  (2),  Ie(Vo)=4zr02eno  (Te/2mn e)  ,/2(/e V0/Te>  is  the  dc 
electron  current,  and  Xo=(Te/4Tze2 no)  is  the  electron  Debye  length.  From  this  expression 
we  may  determine  the  electron  temperature  by  finding  Rac  at  two  or  more  dc  voltages  in 
the  frequency  range  indicated.  Solving  Eq.  (1)  for  Te, 


T.-  e<V‘-V >>  (3) 

In  [RJV2)/R„Ar,)l 

We  also  note  in  passing  that  Eq.  (2)  may  be  inverted  to  solve  for  plasma  potential  given 
an  independent  detennination  of  Te  and  ne.  In  the  work  here  we  have  used  voltage 
differences  from  1  to  3  V. 


b.  Circuit  model 


In  order  to  determine  Te  from  Eq.  (3)  we  must  be  able  to  compute  Rac  from  the  network 
analyzer  measurements  which  provide  Zac(co).  We  consider  the  plasma  as  being  divided 
into  a  quasineutral  core  plus  a  non-neutral  sheath  near  the  probe.  If  the  rf  period  is  much 
greater  than  the  time  for  electrons  to  cross  the  sheath,  the  electrons  remain  in  equilibrium 
with  the  total  voltage,  and  therefore  Rac  — ►  (dlJdV)'1 .  However  there  is  a  capacitance 
associated  with  the  sheath  width  which  lies  in  parallel  with  Rac.  The  circuit  model  we  use 
to  approximate  our  probe  plasma  interface  is  shown  in  Figure  1.  The  primary 
contribution  to  the  total  ac  resistance  arises  from  both  the  real  current  through  Rac  and  the 
displacement  current  through  Cs.  The  net  ac  sheath  impedance  becomes, 


Zac,s(®)  = 


icoC,  +- 


1 


R  +  icoL , 

ac  sh 


(4) 


where  Ls/,  is  the  sheath  inductance  seen  in  Fig  1,  which  we  include  for  completeness  and 
calculate  below.  We  note  that  Cs  will  decrease  as  the  dc  voltage  V0  is  made  more 
negative  and  the  sheath  widens. 
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In  treating  the  region  outside  the  sheath,  we  assume  the  plasma  is  cold  and  collisionless 
and  that  ac  current  is  carried  solely  by  electrons.  The  ac  current  density  then  satisfies 
dJac/dt  =  icoJae  =  e2 Ene/me  based  on  the  electron  momentum  equation.  The  ac  radial 
electric  field  is  thus  given  in  this  region  by  Eac  =  icoJacme/e2ne.  The  ac  voltage  drop  A  Vac 
outside  the  sheath  equals  the  radial  integral  of  Eac  ,  and  since  Jac(r)  =  Iac/4nr 2  in  the 
absence  of  sources  and  sinks,  AVac/Iac  =  icoLp  by  definition.  With  these  assumptions  the 
bulk  plasma  inductance  is  given  by, 


L, 


=  ]dr~ 


m , 


(5) 


4  7tr2e2ne(r) 

Since  the  sheath  region  also  is  characterized  by  a  time-varying  electric  field  induced  by 
the  ac  signal,  the  same  considerations  apply  in  calculating  a  sheath  inductance  or, 

rc 

m, 


E 


'sh 


=  \dr" 


4n :r2e2nes(rY 


(6) 


In  these  expressions  rs  is  the  sheath  radius.  The  sheath  density  profile,  nes(r),  is 
determined  by  solving  the  Poisson  equation  to  determine  rs  and  the  approximate 
sheath/presheath  boundary.  The  plasma  density  ne(r)  in  the  presheath  is  determined  by 
assuming  quasi-neutrality.  With  these  considerations  and  noting  that  the  plasma 
inductance  is  in  series  with  the  sheath  impedance,  the  expression  for  the  total  ac 
impedance  represented  by  the  circuit  of  Fig.  1  becomes, 


Zaci(°)=ZacA(0)  +  i(0LP  = 


P  .  1 
icoC  +  - 


-  +  icoL„ 


R  +  icoL  , 

ac  sh 


or  expanding, 


7  ,  ,  Rac\-(°2LPCS  -l]  +  i(o[Q)2CsLshL  -(Ls  +L  )\ 
Z„c(®)  =  - 


{a>  CsLsh  -1  )-icoCsRa 


After  some  algebra  we  find, 


(7) 


(8) 


R  e(zJ  = 


R. 


(co2LshCs-l)2+(o?RacCs)2 


and, 


(9) 


T  ,  [«>%hLpCs -a>(Lsh  +Lp)}(e>%hCs-l)  +  e>RlCs(a>2LpCs  -1) 

Im(Z  )  = - - - t — - r - r - - -  (10) 

{<o2LshCs-\f  +{<oRacCsf 


We  notice  from  Eq  (9)  that  Re(Zac)  is  a  function  of  co  whereas  Rac  as  seen  in  Eq  (2)  is 
independent  of  co  and  is  a  constant  for  given  plasma  conditions.  Also  Re(Zac)  becomes 
purely  resistive  (=  Rac )  as  co— >  0. 
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c.  Numerical  Algorithm 


Since  the  network  analyzer  returns  Re(Zac)  and  Im(Zac),  and  the  inductances  are 
calculated  based  on  results  of  the  Poisson  equation  solution,  Eqs  (9)  and  (10)  are  two 
equations  in  the  two  unknowns,  Rac  and  Cs. 

In  order  to  detennine  the  inductances  of  Eqs.  (5)  and  (6),  the  Poisson  equation  is 
used  to  construct  the  density  profiles  from  the  spherical  surface  to  pre-sheath. 
Quasineutrality  is  then  invoked  to  detennine  ne(r)  from  the  pre-sheath  edge  to  infinity. 
This  is  accomplished  iteratively  by  specifying  the  sheath  radius  as  a  boundary  condition 
and  integrating  inward  to  the  sphere  bias  voltage  which  will  produce  that  sheath  radius 
for  a  given  plasma  density  and  temperature.  Once  the  bias  is  determined  for  a  given 
sheath  size,  both  the  sheath  and  bulk  plasma  density  profiles  are  known. 

The  algebra  for  solving  Eqs  (9)  and  (10)  is  burdensome  as  they  are  not  easily 
separable;  however,  they  may  be  solved  numerically.  With  the  solutions  for  the  Rac’ s  as  a 
function  of  applied  bias  and  frequency,  Eq.  (3)  may  be  solved  for  Te.  We  show  below 
these  results  and  comparisons  to  independent  measurements  using  conventional  fitting 
methods.  In  the  cases  below  we  show  results  found  by  sweeping  the  spherical  probes 
themselves  both  before  and  after  the  frequency  runs  and  applying  conventional  fitting 
algorithms  which  determine  density,  temperature  and  plasma  potential.6 

III.  Experimental  arrangements  and  technique 

Two  of  the  experiments  described  took  place  in  the  large  volume  portion  of  the 
Space  Physics  experimental  facility  at  the  Naval  Research  Laboratory.  This  section  of 
the  experimental  apparatus  consists  of  a  cylindrical  chamber  of  diameter  2  m  and  length  5 
m.  Argon  plasma  densities  in  this  work  varied  in  the  range  of  107  to  109  cm'3.  Typical 
chamber  pressure  was  1  x  10’4  Torr.  A  typical  electron-neutral  collision  frequency  for 
these  plasmas  is  ven  ~6x  ltfi  s'1  which  is  much  less  than  the  plasma  frequency,  copo  ~6x 
10  s'  (  fPo~  100  MHz).  Neutrals  and  ions  are  at  room  temperature.  The  plasma  is 
created  by  a  recently-redesigned  tungsten  filament  source  biased  to  near  -70  Volts  and 
covering  a  large  portion  of  the  inner  end-plate  surface  area.  A  low-level  axial  magnetic 
field  on  the  order  of  1-2  Gauss  is  provided  by  5  circular  water-cooled  coils  aligned 
axially  in  a  Helmholtz  configuration;  magnetic  field  strengths  near  1  kG  are  available  for 
other  experimental  programs.  Some  of  the  experiments  done  here  were  done  without  an 
external  field.  Electron  density  and  temperature  measurements  are  typically  done  with  a 
cylindrical  Langmuir  probe  which  is  constantly  heated  to  prevent  contamination  buildup5 
or  with  an  emissive  probe.  Further  details  of  the  experimental  configuration  and  the 
general  laboratory  setup  for  different  configurations  are  to  be  found  elsewhere.1314 

The  other  experiment  covered  here  took  place  in  a  smaller  chamber  originally 
designed  and  used  for  dusty  plasma  experimentation.  The  DUsty  PLasma  Experiment 
(DUPLEX)  chamber  is  a  transparent  polycarbonate  cylindrical  vacuum  chamber  which  is 
46  cm  in  diameter  and  61  cm  in  height.16  The  experimental  volume  can  be  accessed 
from  the  top  by  an  aluminum  plate  with  five  vacuum  flanges.  The  cylinder  rests  on  a 
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stainless  steel  chamber  plate  through  which  access  is  provided  to  a  4-in  diffusion  pump- 
based  vacuum  system.  At  the  bottom  of  this  chamber,  separated  by  approximately  50  cm 
from  the  spherical  impedance  probe  used  in  these  studies,  is  a  biased  filament  DC 
discharge  source  which  produced  an  argon  plasma  with  a  typical  density  near  10  cm"  at 
a  pressure  0.1  <  PN  <  0.5  mTorr.  The  biased  discharge  source  voltage  is  in  the  range  of 
about  -50  to  -70  Volts.  The  base  pressure  is  near  .01  mTorr.  The  impedance  probe  is 
cleaned  between  runs  by  heating  with  a  high  power  (~  200  W),  retractable  lamp  which  is 
in  close  proximity  during  the  cleaning  process  and  retracted  during  experimentation. 
Experiments  done  in  this  facility  used  exclusively  the  3  mm  radius  sphere.  The  plasma 
environment  was  characterized  typically  by  higher  electron  temperatures  than  in  the  large 
chamber  environment. 


Three  small  stainless  steel  spheres  of  3  mm,  9.5  mm  and  12.5  mm  radius  were 
connected  to  an  HP8735D  Network  Analyzer  through  50  Q  coaxial  cable  for  the  data 
presented  in  this  work.  The  9.5  and  12.5  mm  spheres  were  mounted  on  a  1/4  inch 
diameter  ceramic  and  steel  support  which  is  connected  to  1/4  inch  diameter  semi-rigid 
copper  50  Ohm  coaxial  cable.  This  improvement  over  previous  versions  of  this 
diagnostic11  was  made  so  that  the  probes  could  withstand  much  higher  temperatures 
while  being  sputter  cleaned.  Primary  cleaning  of  the  3  mm  probe  was  accomplished  with 
a  retractable  heat  lamp  as  described  above. 

For  all  of  the  experiments,  the  determination  of  plasma  impedance  depends  upon 
the  network  analyzer  measurement  of  the  complex  reflection  coefficient,  T(ro).  From  this 
measurement  the  analyzer  returns  as  separate  outputs  Re  Zac(co)  and  Im  Zac(co )  where, 


zac(®)  =  Zo[ 


l  +  r(G>) 
l-r(ru) 


(13) 


and  Zo  (=50  Q)  is  the  internal  impedance  of  the  analyzer.  We  also  note  that  the  ratio  of 
reflected-to-total  power  is  given  by, 

I  |2  P 

rr=—  (14) 

II  p 

1  0 


where  Po  =  Pr  +  Pt  with  PR  and  Pj  the  reflected  and  transmitted  powers,  respectively. 
(The  quantity  1-  |  T  |  2  is  the  nonnalized  transmitted  power  and  this  output  is  also 
available).  The  impedance  from  the  cabling  and  support  is  compensated  through 
instrument  calibration  when  connected  to  a  50  f2  resistor  or  when  calibrated  as  an  open, 
or  short,  circuit.  An  open  circuit  corresponds  then  to  T=l,  a  short  circuit  to  T=-l  and  if 
the  load  impedance  is  a  perfect  match,  T=0.  As  the  change  in  the  complex  reflection 
coefficient  for  the  sphere  in  the  plasma  is  very  small,  this  calibration  is  a  critical  step. 
Care  must  be  taken  to  avoid  unwanted  rf  noise  or  reflections  from  the  chamber  walls  or 
other  nearby  probes.  The  method  is  tested  by  connecting  other  known  resistances  and 
capacitances  to  the  end  of  the  probe  shaft  to  ensure  that  any  error  is  much  smaller  than 
the  changes  in  the  impedance  we  wish  to  measure.  These  considerations  are  covered  in 

1 3 

more  detail  elsewhere. 
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IV.  Experimental  Results  and  Temperature  Measurement  Comparisons 

a.  Applied  bias  levels  and  frequency  range  considerations 


When  applying  Eq.  (3)  to  find  Te  for  an  argon  plasma,  it  can  be  shown  that  biases  could 
range  as  high  as  about  -9  Volts  for  a  1  eV  plasma.  That,  however,  is  not  the  primary 
limiting  factor  in  this  case.  The  bias  levels  applied  to  the  probe  for  the  purpose  of 
calculating  temperature  as  indicated  in  Eq.  (2)  typically  ranged  from  0  to  -3  Volts  in  steps 
of  1  Volt.  The  restriction  on  applicable  bias  levels  is  demonstrated  in  Fig.  2  which 
derives  from  the  Boltzmann  approximation,  which  allows  us  to  express  the  electron 
plasma  frequency  at  the  probe  as, 


^V(V)  =  0),ne 


-AF727I 


pr 


(15) 


where  Te  is  the  electron  temperature,  M  is  ion  mass  and  oopi  is  the  ion  plasma  frequency. 
In  this  expression  AV  is  the  applied  bias  with  respect  to  plasma  potential.  Here,  and 
throughout  the  paper,  we  stress  that  we  only  consider  biases  less  than  the  plasma 
potential.  The  restrictions  both  on  applicable  bias  levels  and  frequency  range  information 
of  Eq.  (15)  can  be  inferred  from  Fig.  2.  These  data  were  taken  with  the  0.95  cm  radius 
probe  in  the  large  vacuum  chamber  where  for  this  run  neo  =  2.3  x  10  cm'  and  the  plasma 
potential  was  +1.9  Volt.  We  use  an  approximate  Te  found  from  conventional  analysis  in 
plotting  Eq.  (15).  The  figure  indicates  decreasing  plasma  frequency  (electron  density)  at 
the  probe  surface  for  increasing  (negative)  bias  as  would  be  expected.  From  the  figure,  it 
is  clear  that  for  a  bias  level  of  AV  =  -3  V,  the  frequency  range  used  to  find  Te  is  restricted 
to  less  than  about  23  MHz.  Since  for  this  argon  plasma,  the  ion  plasma  frequency,^,;,  is 
near  0.5  MHz,  we  are  restricted  to  a  fairly  narrow  range  of  available  frequencies  (i.e., 
( opi«  co  <  cope(r0)).  For  comparison,  the  bulk  electron  plasma  frequency  is  about  136 
MHz.  For  biases  more  negative  than  -3  V,  the  frequency  ranges  become  too  narrow  to 
get  a  sufficient  representation  of  Rac  (which  is  not  a  function  of  co)  and  hence  a  solution 
to  Eq.  (3)  becomes  more  difficult.  Further,  the  curve  of  Fig.  2  serves  to  delineate  the 
cutoff  frequency  for  the  collisionless  resistance  model;  for  a  given  sphere  bias  there  will 
be  no  contribution  of  collisionless  resistance  to  the  net  ac  impedance  for  any  lower 
frequency  than  that  at  the  given  bias  in  the  figure.  At  higher  frequencies  the  net 
impedance  will  contain  contributions  from  the  collsionless  resistance  as  well  as  from  the 
lower  frequency  impedance  treated  in  this  paper. 

b.  Data  and  determination  of  Te 

The  first  data  set  we  analyze  is  based  on  the  3  mm  radius  stainless  steel  sphere  in  the 
DUPLEX  chamber  described  above.  We  demonstrate  in  graphical  form  the  construction 
method  for  determining  Te  from  the  data  taken  by  the  network  analyzer.  For  a  portion  of 
the  later  data  using  different  spheres  of  varying  diameter  in  a  larger  chamber,  we  show 
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only  the  results  of  the  measurements  along  with  comparisons  to  determination  of  Te  by 
standard  fitting  routines  used  with  the  spheres  operated  as  conventional  Langmuir  probes. 

Plotted  in  Fig.  3  is  the  comparable  plot  to  Fig.  2  for  the  3mm  probe  data.  From 
comparison  of  the  two  plots  one  can  see  the  strong  impact  of  Te  on  the  allowed 
frequencies  before  cut-off.  i.e.,  whereas  in  Fig.  2  for  Te  ~  1  eV,  a  voltage  bias  of  -4  Volts 
places  the  cutoff  at  near  14  MHz,  in  the  case  of  Fig.  3  at  a  comparable  density  but  at  Te  ~ 
3  eV,  it  is  not  possible  to  analyze  the  data  to  determine  Te  at  this  bias.  Even  for  a  bias  of 
-3  Volts  the  frequency  range  allowable  is  severely  restricted  i.e.,fpi  ~  0.6  MHz  and  fpe(ro) 
~1.8  MHz.  Therefore  for  this  temperature  and  density  we  are  restricted  to  biases  no 
larger  than  -  -3  Volts. 

Fig.  4  is  a  plot  of  Re(Zac)  versus  frequency  for  probe  biases  varying  from  0  to  -4  volts  for 
the  3  mm  radius  sphere.  This  example  shows  the  frequency  dependence  of  the  function  as 
seen  in  Eq.  (9).  We  contrast  this  to  Fig.  5  where  the  results  for  Rac  vs  applied  frequency 
are  plotted.  Rac  is  averaged  over  the  frequency  ranges  and  this  mean  value,  for  a  given 
applied  bias,  is  plotted  as  a  straight  line  and  used  to  determine  Te  using  Eq.  (3).  Fig.  6 
shows  Te  versus  the  voltage  range  bin  chosen,  along  with  the  mean  value  (black  solid 
line)  and  the  value  of  Te  determined  by  our  Langmuir  probe  fitting  routines  (red  solid 
line).  As  an  example  if  we  select  voltage  bins  which  include  voltages  varying  from  0  to  - 
3  Volts  in  steps  of  1  V,  this  allows  6  different  calculations  of  Te  corresponding  to  the 
voltage  bins:  (1)  0,-1  V  ;(2)  0,-2  V;  (3)  0,-3V;  (4)  -1,-2V;  (5)  -1,-3V;  and  (6)  -2,-3V.  In 
the  case  of  Fig.  5  we  have  included  biases  from  0  to  -3  Volts  even  though  the  final 
voltage  bin  6  (-2,-3  V)  is  again  restricted  in  the  number  of  usable  frequencies  according 
to  Fig  3. 

The  results  for  the  9.5  mm  probe  in  the  large  plasma  chamber  environment  described 
above  are  seen  in  Figs.  7  and  8  for  B=0  and  in  Figs  9  and  10  for  B=1  Gauss.  Figs  1 1  and 
12  show  only  Te  for  the  12.5  mm  radius  probe  under  conditions  differing  only  slightly  in 
Te  and  ne 

As  a  final  example  we  show  a  data  set  taken  for  the  12.5  mm  probe  at  low  density  (ne  ~ 
1.2  x  107  cm-3  )  and  higher  temperature  (Te  -  2.6  eV)  using  an  He  plasma.  The  plasma 
potential  is  near  3.8  Volts.  Because  of  the  low  electron  density  we  are  again  restricted  in 
terms  of  the  frequency  range  in  order  not  to  include  any  contribution  from  the  resonant 
CR  as  covered  above.  The  biases  vary  only  from  0  to  -2  V.  The  cutoff  frequency  is  near 
10  MHz  but  we  further  restrict  that  frequency  to  6  MHz  to  be  safely  below  this  range. 
Instead  of  using  a  value  of  Rac  which  is  averaged  over  the  available  frequency  ranges  as 
in  the  earlier  plots,  here  we  plot  Te  vs  applied  frequency  for  each  of  the  bias  voltage 
ranges  separately  in  Fig.  13.  The  plot  shows  3  values  for  each  frequency  applied 
corresponding  to  the  voltage  bins  for  a  0  to  -2  Volt  range.  As  in  the  data  above  the 
standard  deviation  is  taken  for  the  points  about  the  mean  value  of  calculated  (theoretical) 
temperature. 


V.  Summary 


We  have  introduced  a  method  of  determining  electron  temperature  for  small 
spherical  probes  using  results  of  a  network  analyzer  determination  of  impedance.  The 
method  does  not  depend  on  fitting  routines  and  hence  does  not  require  a  specific  form  to 
be  used  while  interpreting  an  I-V  characteristic.  However,  it  does  require  estimates  for 
Lsh,  Lp  and  Cs.  The  method  is  robust  in  the  sense  that  it  not  affected  by  probe 
contamination  and  does  not  formally  depend  upon  the  assumption  of  a  single  Maxwellian 
distribution.  Finally,  the  number  of  independent  determinations  of  Te  can  be  increased  by 
using  finer  resolution  voltage  bins  than  we  have  used  in  this  work. 

As  can  be  seen  in  the  comparisons  above,  temperatures  obtained  by  the  method 
tend  to  differ,  in  some  cases  more  and  some  less,  from  those  determined  by  a  fit  using 
conventional  analysis.  Generally  the  temperatures  in  the  large  chamber  using  the  new 
method  are  higher  than  those  produced  by  the  conventional  fitting  procedures.  We  plot 
in  Fig.  14  the  values  for  Te  determined  by  the  different  methods  in  this  paper  versus  each 
other.  Although  it  is  clear  that  generally  the  temperatures  measured  by  the  impedance 
method  outlined  here  are  higher,  the  agreement  appears  reasonable  between  the  two 
methods.  One  possible  source  of  error  in  conventional  I-V  analysis  involves  the  ion 
contribution  to  the  total  current,  particularly  for  negative  bias  levels.  For  example,  the 
existence  of  a  small  component  of  thermalized  hot  electrons  can  have  a  major  effect  on 
the  determination  of  ion  current  in  the  ion  saturation  region.  This  holds  true  if  a  two- 
temperature  distribution  is  present  and  not  accounted  for.  If,  due  to  the  presence  of  the 
small  hot  electron  component,  the  ion  saturation  current  is  estimated  to  be  smaller  than  it 
actually  is,  this  can  influence  the  calculation  of  the  bulk  plasma  temperature  depending 
on  the  particular  analysis  method  employed.  In  a  recent  work17,  we  show  that  a 
surprisingly  small  component  of  higher  temperature  can  have  a  dramatic  effect  on  ion 
saturation  current  determination. 
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VI.  Figure  Captions 


Figure  (1)  -  Circuit  schematic  analogy  to  probe/sheath/plasma  used  to  estimate 
impedance. 

Figure  (2)  -  Plasma  frequency  at  the  probe  surface, /L/ra)  ,  vs  applied  bias  voltage,  Vbias, 
wrt  to  plasma  potential,  Vpiasma.  Ne  =  2.3  x  10'  cm'  ,  r0  =  .95  cm 

Figure  (3)  -  Plasma  frequency  at  the  probe  surface,  fpe(ro)  ,  vs  applied  bias  voltage,  Vbias, 

o  •* 

wrt  to  plasma  potential,  Vp/asma.  Ne=  5.3  x  10  cm'  ,  ro  =  .33  cm 

Figure  (4)  -  Re(Zac)  vs  frequency  for  5  different  bias  voltages.  Ne  =  5.3  x  10  cm'  ,  r0  = 
.33  cm 

Figure  (5)  -  Rac  vs  frequency  and  averages  for  4  different  biases.  (V=0,  black;V=-l, 
red;V=-2;green;  V=-3  blue)  Ne  =  5.3  x  108  cm'3 ,  ro  =  .33  cm 

Figure  (6)  -  Te  fit  (red)  and  Te  average  (black)  vs  bias  voltage  range  bin:  (1)  0,-1  V  (2)  0,- 
2  V,  (3)  0,-3  V  (4)  -1,  -2  V  (5)  -1,-3  V  (6)  -2,-3  V.  Te(average)(black)  =  2.63  eV,  a  = 
0.426  eV  Ne  =  5.3  x  108  cm3  r0  =  .33  cm. 

Figure  (7)  -  Rac  vs  frequency  and  averages  for  4  different  biases.  (V=0,  black;V=-l, 
red;V=-2;green;  V=-3  blue)  .  Ne  =  2.3  x  108  cm'3 ,  r0  =  .95  cm 

Figure  (8)  -  Te  fit  (red)  and  Te  average  (black)  vs  bias  voltage  range  bin:  (1)  -1,-2  V  (2)  - 
1,-3  V,  (3)  -1,-4  V  (4)  -2,-3  V  (5)  -2,-4  V  (6)  -3,-4  V.  Te(average)  =  0.73  eV,  a  =  0.11 
eV,  Ne  =  2.3  x  108  cm3  r0  =  .95  cm. 

Figure  (9)  -  Rac  vs  frequency  and  averages  for  4  different  biases.  (V=0,  black;V=-l, 
red;V=-2;green;  V=-3  blue)  .  Ne  =  3.38  x  10s  cm'3 ,  ro  =  .95  cm,  B  =  1  Gauss 

Figure  ( 1 0)  -  Te  fit  (red)  and  Te  average  (black)  vs  bias  voltage  range  bin:  (1)  -1,-2  V  (2)  - 
1,-3  V,  (3)  -1,-4  V  (4)  -2,-3  V  (5)  -2,-4  V  (6)  -3,-4  V.  Te(average)  =  1.06  eV,  a  =  0.17 
eV,  Ne  =  3.38  x  108  cm3  r0  =  .95  cm.  B  =  1  Gauss 

Figure  (11)  -  Te  fit  (red)  and  Te  average  (black)  vs  bias  voltage  range  bin:  (1)  0,-1  V  (2) 
0,-2  V,  (3)  0,-3  V  (4)  -1,  -2  V  (5)  -1,-3  V  (6)  -2,-3  V.  Te(average)  =  1.85  eV,  a  =  0.076 
eV  Ne  =  5.2  x  107  cm3  ro  =  1.25  cm. 

Figure  (12)  -  Te  fit  (red)  and  Te  average  (black)  vs  bias  voltage  range  bin:  (1)  0,-1  V  (2) 
0,-2  V,  (3)  0,-3  V  (4)  -1,  -2  V  (5)  -1,-3  V  (6)  -2,-3  V.  Te(average)  =  1.59  eV,  a  =  0.12 
eV  Ne=  3.6  x  107  cm3  r0  =  1.25  cm. 

Figure  (13)  -  7)  lit  (red)  and  Te  average  (black)  vs  frequency  and  averages  for  5  different 
biases,  a  =  .18  eV,  Te(average)  =  2.59  eV,  Ne  =  1.27  x  107  cm'3 ,  r0  =  1.25  cm,  B  =  1 
Gauss 
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Figure  (14)  -  Te  fit  vs  Te  (impedance)  for  6  different  experimental  series. 
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